Get default configuration
FUNCTION get_default_config() RESULT(config) TYPE(config_type) :: config ! Default values are already set in the type definition ! No need to explicitly set them as they're defined in the type config = config_type(pivot_tolerance=TOL_PIVOT, & convergence_tolerance=TOL_CONVERGENCE, & residual_tolerance=epsi, & max_iterations=kmax, & block_size=64, & use_openmp=.TRUE., & use_blas=.TRUE., & preallocate_workspace=.TRUE., & workspace_size=1000, & enable_debug=.FALSE., & enable_timing=.FALSE., & log_file="nafpack.log", & default_direct_method="A_LU", & default_iterative_method="Gauss_Seidel", & default_preconditioner="ILU") END FUNCTION get_default_config